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ABSTRACT 


Discretization  formulas  for  the  digital  simulation  of  the  dif- 
fusional  processes  associated  with  a  variety  of  electrochemical  processes 
at  a  spherical  electrode  are  given.  The  formulas  are  based  on  orthogonal 
collocation  techniques.  Several  mechanisms  are  represented,  and  the 
technique  is  demonstrated  for  chronoamperometry ,  chronopotentiometry , 
and  cyclic  voltammetry. 
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I 

Orthogonal  collocation  techniques  for  the  simulation  of  second 
order  partial  differential  equations  have  been  demonstrated  for  a 
variety  of  problems  in  electrochemistry  [1-7],  chemical  engineering 
[8-11],  differential  scanning  calorimetry  [12],  and  other  fields  [13]. 
The  advantages  over  other  methods  have  been  described,  but  noteworthy 
is  generiilly  increased  accuracy  for  decreased  computational  effort 
[1-7]  . 

In  electrochemical  diffusion  problems,  the  algorithms  needed  to 
compute  concentration  profiles,  chronopotentiograms,  chronoamperometric 
responses,  spectrophotometric  responses  of  individual  species,  cyclic 
voltammograms ,  and  chronocoulometric  responses,  all  to  planar  electrode 
configurations,  have  been  given.  This  paper  deals  with  some  of  the 
same  experiments  applied  to  a  common  configuration,  the  static 
spherical  electrode  (Figure  1) . 

The  programs  for  solving  the  equations  described  herein  have 
been  given  [7].  Basically,  the  orthogonal  collocation  technique  is 
implemented  by  the  following  steps: 

1.  Make  the  Fick's  laws  equations  suitably  dimensionless. 

2.  Discretize  the  resulting  equations  at  the  roots  (zeros)  of  an 
orthogonal  polynomial  of  suitable  symmetry.  The  resulting  set  of 
equations  are  now  first  order  ordinary  differential  equations. 

3.  Integrate  the  set  of  simultaneous  differential  equations  to  obtain 
the  set  of  concentration  profiles  with  respect  to  distance  from 
the  electrode  surface  and  with  time. 

4.  Use  these  concentrations  in  the  suitably  discretized  equations 
to  obtain  the  desired  electrochemical  parameter  (current,  absor¬ 
bance,  etc.). 
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Chronoamperometry 

We  will  only  consider  diffusion  on  the  solution  side  of  the 
electrode-solution  interface.  For  a  chronoamperometric  experiment  on 
the  simple  charge  transfer  mechanism  at  a  spho’-’cal  electrode 

ne~ 

A  ~ ^  B  (1) 

we  have  the  following  imposed  conditions: 


fAlR,0  ~  {AJco,T  lA°J 

[A]  m  =  [B]  =  [B]_  n  =  0 

rQ,I  oo,T  R,U 

(2) 

/<5  f  A ] \  / 6  f  B]\ 

\iR  /R=r0  "  T?R  7R-r0 

where  R  is  the  distance  parameter  (radial)  from  the  center  of  the 
spherical  electrode  and  rQ  =  the  radius  of  the  electrode.  Since 
interpolation  methods  using  orthogonal  polynomials  are  simpler  to  use 
in  a  [0,1]  interval,  we  will  define  M  as  some  distance  in  the  R 
direction  such  that  no  diffusion  effects  are  experienced  there  during 
the  time  frame  of  the  experiment.  The  differential  equation  for  the 
A  species  is: 

xr*i  D.62[A]  2D.  . . . , 

6 1 A  ]  _  A  _  ,  A  5  [  A  ]  #  \ 

6T  £r2  '  R  6R "  1 

For  dimensionless  variables,  we  choose 


t 


DT 


(N-r0) 


r 


R-r 

M~-r 


0 

0 


(4) 


Wc  note  then  that 
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dT  = 


(M-r0) 


dt,  and 


dR  =  (M-rjj)dr 

Substitution  of  the  new  variables  (4)  into  (3)  gives 

„  DA|A<,1s2cA  +  (  2DA  /^i^A 

(M-rQ)26t  (M-rQ)26r2  (M-rQ)r+r0  (M-rQ)6r 

which,  after  simplification,  leads  to 

6c.  62c.  2  6c 

_ A  _  _ A  _  _ A 

6t  6r2  r  +  g  6r 

where 


(5) 


(6) 


(7) 


The  general  spatial  derivatives  in  terms  of  their  polynomial  inter¬ 
polation  coefficients  [1]  are  given  by 


N+2 

I 

j-1 


Q(X.  t) 

J  » 


N+2 

=  E  B.  .  Q (X .  t) 
j-1  1,3  3' 


l 

where  N  is  the  order  of  the  approximation  polynomial  chosen, 
represents  the  roots  of  that  polynomial  collocation  points, 
j  and  are  given  by  [10] 


(8a) 

(8b) 

the  Xi 
and  the 


4 


•V 


dX^/dX 

dX°/dX 


dxJ+2/d>: 


dxJ/dX 


dX^+2/dX 


dx^2/dx 


-1 


(  9) 


and 


|Bij>  * 


-d2xJ/dX2  d2xJ/dX2 

d2>:2/dx2 


d2X°+2/dX2 


d2X^+2/dX 


.2  N+2  ,  ,Y 
d  XN+2^dX 


-1 


(10) 


The  matrix  [S]  is  given  by 


[S]  = 


rY0  Y1 
X0  X0 

X1  ■' 
X° 

*-XN+2 


-  X 


N+2  n 
0 


.N+2 

N+2 


(ID 


The  algorithms  for  finding  the  roots  of  any  orthogonal  polynomial  and 
the  corresponding  [A^]  and  (B^]  matrices  have  been  given  previously 
19].  The  N+2  terms  arise  because  X=0  and  X=1  are  included  as  roots  to 
the  orthogonal  polynomial.  The  concentration  profile  will  be  eval¬ 
uated  at  the  collocation  point,  i.e.  at  coencentric  spheres  surrounding 
the  spherical  electrode.  Each  sphere  has  radius  r=X^ ,  X^ »  X^,  etc., 
i.e.  at  the  collocation  points. 

Equation  (7)  is  thus  written  in  terms  of  (8a)  and  (8b)  at  the 
discrete  points  X^: 


dc 


A 
d  t 


N+2 


Nt  2 


^  Z  13.  . c  (r .  1)  +  — Z  A.  .c  (r  . 
ri-Xi  j  =  l  1J  A  3'  ri+l'  j=-l  1J  A  3 


t) 


(12) 
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Thus  we  have  N+2  simultaneous  first  order  differential  equations  to 

solve  for  the  N+2  c.(r,  t)  unknowns.  Any  one  of  several  methods  may 

A  3 , 

be  used  for  solving  this  system  of  equations  {1,4,9], 

Inserting  the  boundary  conditions  { 2  )  further  defines  the 
electrochemical  method  used,  and  reduces  the  number  of  equations  from 
N+2  to  N.  This  is  done  by  first  expanding  equation  (12)  partially, 
exposing  the  boundary  points  at  r^=X^=0  and  rN+2=XN+2=* : 

dCA 
dt 


=  BiriCA(Oft)  +  Bi/N+2cA(l,t)  +  z  B.jcA(rj/t)  + 
ri  x  D~ £ 


ri*e  [AincA 


(0,t)  ♦  AlfH+2cA(l,t)  *  ^  Ajj 


We  have  from  the  original  boundary  conditions  (  2  )  thcit  cA(0,t)  =  0 


and  cA(l,t)  =  1,  so  that  (14)  becomes 


r  .=X . 
x  1 


N+i.  r  n+i  "j 

=  B.  „.-+  I  B .  .c,  (r  .  t)+ — — r  A .  -+  Z  A.  .c-  (r  .  t) 

l,N+2  j=2  13  A  3,  ri+E[  i.N+2  _.  =  2  13  A  3,  J 


So  we  have  one  equation  to  solve  at  each  of  the  N  interior  collocation 
points  ,  i  =  2,  .  .  .N+l  for  the  N  unknown  cA(r^  t). 

The  concentration  profile  for  the  B  species  is  derived  similarly, 
and  leads  to  identically  the  same  result  as  equation  (12): 


am  r 

i  i  J  J 

Parti.il  expansion  leads  again  to  equation  (13)  with  c^tr^  t)  substituted 

for  the  c. (r.  t) .  However,  even  though  it  is  true  that  one  boundary 
A  3 » 
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condition,  Cg(i,t)  =  1  may  be  immediately  substituted,  the  other  boundary 
condition  Cg(0,t)  is  not  known  explicitly.  This  problem  is  easily 
resolved  in  collocation  methods.  We  observe  that  the  flux  relation 


(16 


We  will  let  D^=D^  here  for  convenience  even  though  it  presents  no 
problem  in  setting  up  the  equations  if  we  do  not  [1,4].  The  flux 
relation  is  replaced  by  equation  (8a)  on  each  side: 


N+2 

1  A1  iCA(ri  = 

i=l  -*■  •  J  A  D  * 


N+2 

-  Z  A  c  (r  t) 
i=l  1 •  J  B  3  ' 


(17 


Partially  expanding  so  that  known  conditions  may  be  inserted,  wc  have 


N+l 


A1,1CA(0' t)+Al,N+2CA(1't)+i^2  Al, jCA(rj,fc) 


-[Al,lCB(rj,t)+Al,N+2CB(1't)+J2  Al,jCB(rj,t)] 


N+l 


(18) 


Inserting  known  quantities,  we  deduce 


N  1 1 


N+l 


AUN+2*i;:i  t  (19) 


Solving  for  our  unknown  boundary  condition. 


N+) 

c  (0,t)  -  U  +  T.  C  (c  (r  t)  +  c  (r  .  t)) 
u  ^.*2  1 9  .i  **  J  •  &  J  * 


(20) 


where 


7 


Al,N+2 

U  =  - -  ,  and 

*1.1 


A,  . 

C  = 

C10  A 


We  can  now  substitute  this  explicit  value  for  Cg(0,t)  in  equation 


(15),  along  with  cD(l,t)  =  0  to  obtain  the  equation  necessary  to 

O 


represent  the  concentration  profiles  for  B: 


dc 


B 


dt 


N+l  N+l 

=  V,  +  Z  D . . c. (r  .  t)  +  Z  E..c0(r.  t) 


r.=X. 
j  l 


i  j_2  ID  A  j, 


j=2 


ij  B  j 


where 


2A.  .U 

V.  =  B.  .T  +  — 
l  i,l  r.+6 


D .  .  =  C .  . 

ID 


B. 


2A .  , 
i ,  1 


i,l  ri+g 


E.  .  =  B.  .  +  C.  . 
ID  1D  i'D 


G. 


2A.  . 
l ,  1 


i , 1  ri+B 


2A .  . 

+  — il 


r  .+£ 

l 


Again  we  have  N  equations  in  the  N  unknown  c  (r.  t)  to  solve 

B  3 , 


simultaneously. 

Table  1  compares  results  obtained  for  the  current  obtained  by 
equation  (33)  and  the  analytic  solution  for  this  case,  which  is 


nFADA1/2lA°)  nFAD^lA0] 

i  = - ; — 7-r. — + 


We  use,  by  example,  =  1.  >:  10  ’’  cm2s  ^ ,  (A0)  -  10  ®  moles  cm 


2  4 

A  =  nr^,  Tq  -  0.1  cm,  and  n  *  1.  £  for  the  example  was  10 


For  the  mechanism, 


<21)| 


(22 


■J" 


8 


I 

I 


i 

i 

i 


.  ! 


I 


< 

J 


;  ( 


I  ! 

8 


4 

I  1 


A  ;=2  B 


k 

— >  c 


the  profile  for  the  A  species  is  the  same  as  equation  (14)  .  The  B 
species  has  a  different  profile  due  to  loss  of  material  through  the 
chemical  reaction,  i.e. 


6[Dj 

6t 


+ 


2n  MB] 

— - k  [B] 

R  6R 


(23] 


After  discretization  to  the  collocation  points  and  insertion  of  the 
dimensionless  variables  (  4  ) ,  we  have 


dc 


B 


dt 


N+2 

=  I  B 
r . =X .  j=l 

l  l  J 


2  N+2 

i.jVq.*1  +  VB"  jh  Ai.j°B,r3-t)  '  “CB(ri,U 


(24: 


where 

a 


k  (K-r 


0 


2 


Continuing  precisely  as  for  the  simple  electron  transfer  case, 
except  with  the  addition  of  this  last  chemical  kinetic  term,  we  find 
that 


N+l 


N+l 


V,;,  DijcA<rj,t)*iL  EijcB(rj,t»-acBlrj,t) 

3  z  j  z 


For  the  catalytic  mechanism 
k 

A  7—  B  — »  C  I  A 


(25 


(26 


both  the  A  and  B  profiles  are  modified  by  the  presence  of  the  chemical 
reaction.  The  differential  equations  are 


MA]  _  <S2(A]  2DA  6  r A 1  .  ,  , 
""fit-  ■  da~~2"  +  ir nC  4  K|B] 

o 

MB]  _  n  <52  [  B]  2DA  MB]  _  .  .  , 
"6“  “  DB~“2“'  +  K  '  ok" 

OK 


(27 


The  boundary  condi tions  arc  exactly  as  the  previous  cases,  so  that 
discretization  is  the  same.  We  find  that 


] 


mrrammm 


dc. 


dt 


N+l 


=B 


N+l 

+  E  A.-c^tr.  t)]+ 


ri=Xi 


,N+2+.i2  BijCA{rj,t)+FT+3tAi/N+2+.?2AijCA<rj 


aCB(rj,t} 


dc. 


dt 


N+l 


N+l 


r  .=>: 

i 


+  E  D .  .c  (r  .  t)+  E  E..c  (r.  t)-ac  (r.  t) 
j„2  A  3>  j=2  b  B  3' 


again  letting  :  =  D_ . 

J  j\  B 

To  calculate  the  current  at  any  time  for  the  mechanisms  listed 
above  under  the  chronoamporomet ric  experiment,  we  note  the  following: 


H  dR  /r=0 


i  =  nFAD 


where  A  is  the  electrode  surface  area  and  F  is  the  Faraday  constant. 
The  flux  tern:  is  convetted  to  dimensionless  terms  by  insertion  of 
the  newly  defined  variables  (4) . 


i  = 


nFAD 

M 


D[A°]  (6ZA 

~r0  \dr  /r1=X1=0 


(28) 


(29) 


(30) 


(31! 


The  flux  term  may  then  be  discretized  to 


N+l 


.=nFA2J^i  [AlflcA(0,t)+A1(B+JcAO,t)+  t 
U  J 


(321 


or 


.  nTAD | A° ]  r.  .  Nt1  .  „  .  .  . , 

(Al,N+2  +  .T=2  Al,jCA(lj,t)1 


(3: 


after  insertion  of  the  known  values  of  c.(0,t)  and  c  (l,t). 

il  s\ 

So  the  current  i  may  Lo  calculated  at  any  time 

(M~ro) 

T  -  - — -  by  insertion  of  the  calculated  concentrations  from 


equations  (14)  or  (28)  into  equation  (33)  .  If  there  are  ad¬ 
ditional  electroactive  species  in  the  mechanism,  their  contribution 
is  simply  added  into  equation  (33),  since  generally, 


i 


nFAD [ A° ] 


M-r 


0 


m 

Z 

n=l 


0 


where  m  is  the  total  number  of  electroactive  species. 

The  charge  or  spectroelectrochemical  absorbance  can  be  cal¬ 
culated  by  use  of  Gaussian  quadrature  weights  [4,9,10]  for  the 
integration  of  the  current  in  time  or  the  concentration  profiles  in 
time . 


Cyc]  ic  Volt  amine  trie  Response 


Several  modifications  have  to  be  made  to  the  method  of  Speiser 
and  Rieker  (6j  to  develop  the  cyclic  voltanunetric  responses  for  the 
spherical  electrode.  We  define  new  dimensionless  time  and  distance 
variables 


and 


t  = 


r  = 


DTa 


(M-rQ) 


(R-r0)a 


1/2 


M-r, 


where 

a  =  nFv/R'T' 

with  v  =  the  potential  sweep  rate  in  volts  s  R'  the  gas  constant, 


(M-r  )dr 
aR=  0_ 

a 

For  the  simple  electron 

transfer  reaction 

(  1  ) ,  the  differential 

equation  ( 6  )  becomes 

aD[A° ] 6c.  aD [A° ] 62c. 

2D [ A° ] 

(M-r  )“6t  (M-rQ)  6r  ( (M-rQ) r+ rQ) (M-r0)  fir 


which  simplifier,  to 


6c  , 

,2  0  1/2 

i>  c.  2a 

6c. 

A 

A  + 

A 

fit 

2  + 

6r  r+6 

6r 

The  derivation  now  proceeds  as  described  by  Speiser  and  Rieker 


The  cyclic  voltamme trie  boundary  conditions  are 


lAlR;0  "  [Al-,t  =  lA°] 


[B]R;0  "  [Bl~,t  "  ° 


=  MM 

R=0  ^  R=0 


TgT^-;=exp[nF/R,l"(E®/B  -  Ell  -  0A/BSjtt) 
U  /  M 


Inserting  the  potential  sweep  rate  (14)  v,  we  have 
®A/B  =  t’XpfniyR'T -  Estakt)J 


Sa(T)  - 


cxp(-nPvT/R’T' )  « cxp(-aT)  for  T  1  T^ 
exp  (at  -  2aT.)  for  T  T  .>  2T 
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where  T  is  the  time  when  the  forward  potential  sweep  is  reversed. 

We  convert  equations  (39  )  to  dimensionless  form  normal¬ 

ized  to  t  and  M: 


cA(r,0)  =  cA(l,t)  =  1 
cB (r , 0)  =  cB  (1 , t)  =  0 


(41) 


cA(0,t) 

^(oTtT  "  eA/Bsx(t) 


Discretizing  using  equations  (8a)  and  (8b)  on  (38) ,  we  have 


N+2 

=  I  B .  .c  (r .  t)  + 


2a1/2  N+2 


^  -jf,  BijcA(rj,t)  +  vr 


or  after  partial  expansion 


uw.  x 

=B.  1  c  (0  / 1)  +B .  (l,t)+  Z  B .  .c  (r .  t)  + 

at  _y,  a,l  A  i  ,N+2  A  .  »  1}  A  j , 

ri  i  3  ,t 

2a1/2  N+l 

rT+fT  lAi,lCA(0,t)+Ai,N+2CA  (1't)+^2AijCA(rj  ,t}  1 


Inserting  the  explicit  boundary  conditions,  we  have 


6T  T  CO.t) ♦** ^ A jCA«r J  t)  +  |S_rIAii!cA(0,U 

r.-xj  3-2  i 


N+l 
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Identical  treatment  for  the  B  species  yields 


dc 


B 


dt 


N  +  l 


=B  ,c  (0,t)  +  i  B..c  (r.  t)  +  p 
r  .-X .  1,1  B  j--2  13  B  3'  ri 


1/2 


[A,  .c  (0,t)  + 


+0  1  i  ,  i  B 


l  l 

N+l 

l  AiicB(ri  t)5 
j=2  13  B  3' 


(45) 


We  discretize  the  flux  relation  in  (41)  to  get  an  equation  relating 
c^(0,t)  and  c  (0,t),  presently  implicit  boundary  conditions. 


N+2 


N+2 

-  E  A.  •  c  (r  .  t) 


j-1 


ir)  E'  3 


(46)1 


Rieker  and  Speiser  [6]  showed  that  after  expansion  of  equation 
(46),  one  may  then  substitute  in  equation  (44)  and  arrive  at  the 
following  result  giving  the  time  (potential)  dependence  of  cA(0,t) 


and  cB (0, t) : 


N+l 


0a/rS\ 

cA(0,t)  -  -  A1';iTl+eA/BSA  (t)  )  Al,N+2+j^2Al,j(GA(rj,t)+CB(rj,t)) 

.  N+l 

CB(0,t)  =  ”  Alfl[l  +  C'A/BSA  (t)  )  Al,N+2+^2Al,  j  (CA(rj,t)+CB(rj,t)  } 


Inserting  these  time  dependent  boundary  conditions  into  the  dis¬ 
cretized  equations  (44)  and  (45)  yield  the  2N  simultaneous  first  order 
differential  equations  necessary  to  solve  for  the  conccntJ ation  profiles 
in  time: 


(47a 


dt 


N+l  N  +  l 

«W.  +  E  F . . c . { r .  t)  +  E  G..cn(r.  t) 
r.=x.  1  j=2  13  A  3 '  j=2  33  B  3' 
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(48) 


dcl 

dt 


N+l 


N+l 


=  X  +  Z  H  c  (r  >  +  Z  l..c  (r  > 
r^=Xi  j=2  j=2  13  B  D*t 


with 


w  _  _  eA/BSl  (t) _ .  +2ayi  , 

i  A,  ,  [l+eA  /BS,  (t)  ]  1  i,lAl,N+2  r  h  B  Ai,lAl,N+2J 


l,ll  uA/B^A 

1/2 


B. 


2a‘ 


A. 


i*N+2  ri+iJ  i,N+2 


F.  .  =  B.  .+- 


2a 


1/2 


A.  .  - 


1/2 

23  Ai, 1A1, jgA/BSX (t) 


ID  ID  ri+P  ij  (V6)Al,lf1+0A/BSX(t)] 


Bi,lAl,-jeA/BSA(t> 

Al,lll-eA/BSX(t)1 

c  _VBSA(t) 

Gij  '  '  “ 


AlflU+6A/BSX(t)1  fBi'lAl'D+(ri+B) 

1/2 


2a1/2A.  .A  . 


Wi  A,  n [l+6a/nS1 (t) ] [Bi,lAl,N+2+ri+e  Ai,lAl,N+2 


H 


ID 


*1,1 1  A/BA 
Al,ltl+eA/BSX(t)1 tBi,lAl,j  '  ri+B~ 


2a1/2A.  A,  . 
1*1  1*D 


2  a 


1/2 


2a1//2A.  ,A.  . 
1*1  1  *  D 


"  Bij  +  TlTT  AiD  '  (Vft,Alflll+eA/BSX(t)1 


B.  ,  A.  . 
1*1  1 »  D 


Al,l(1+0A/ESX(t)] 


The  2N  equations  (48)  are  solved  simultaneously  as  before 

to  find  the  cR  (r  .  t  )  and  c  (r.  t) • 

A  D  *  B  3  * 

The  current  response  is  represented  similar  to  the  chronoamp- 
erometric  cast: 


(49) 


or  in  the  dimensionless  space 

1/2 


1  = 


nFADA[A 


(M 


A°]a^  /dcA\ 

-r0)  \dr  yr^=x^=0 


Discretizing  the  flux  and  inserting  the  boundary  conditions  leads  to  (6) 


dc 


A 


dr 


N+l 


+  £  A  • c  (r .  t)-0*  /DS,  (t) 


N+l 


1  A1  iCB(r-i  t)] 
j_2  1 > J  B  3  » 


So  that 


l  = 


nFADA[A°]a 


1/2 


N+l 


w=V 


.[ _ + 


^  -;cji  ^r-;  /0S,  (t) 


1  +  0A/BSX(t)  l<N+2  j,2  Id  A  j-  uA/B°X 


N+l 


The  concentrations  c  (r.  t)  and  c  (r.  t)  found  at  each  time 

A  1  t  B  X  t 

2 

T  =  (M-rn)  t/aD  are  used  in  equation  (52)  to  determine  the  current 


as  a  function  of  time  and  hence  potential  by  the  relation 


E  = 


(M-rQ)  tR  'T ' 


nFD 


ESTART  +  2vtx 


{M-rQ) ‘tR'T’ 
nFD 


t  £  t 


t  i  t  i  2t 


where  t  is  the  dimensionless  potential  sweep  reversal  time 
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Chronopotentiomctry 

The  chronopotentiometric  case  for  planar  electrodes  was  described  in 
a  previous  paper.  For  a  chronopotentiometric  response  to  the  simple 
reversible  charge  transfer  reaction,  we  must  consider  the  following 
conditions : 


II 

o 

% 

u 

< 

[B]r,0  = 

'b1»,t  =  0 

E  +  ea/b  * 

■BIT*.  ln[B] 

nF  r0 

/fi [A  ]  \ 

/firB]\ 

V  6R  Jr 

V  6H  )r 

v  /ro 

/ro 

i  =  nFAD 

^  =  Constant 

\  6r  /R'r0 

The  new  dimensionless  parameters  are 


t  =  T/t 


r 


R-r 

M-r 


0 

0 


with  t  the  transition  time  of  the  experiment.  Insertion  of  these 
values  into  the  dif ferential  equation  (3)  and  simplification  yields 


2 

6c.  6  c,  ~  fie. 

_ A  _  _ A  ,  2a  _ A 

fit  a  fir  r+3  fir 


where 


a 


Dr _ 

(M-r0)2 


The  conditions  (54-55)  in  dimensionless  form  become 
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cA(r,0)  =  cA(l,t)  =  1 


cB(r,0)  =  cfi ( 1 ,  t )  =  0 


P),  ■  PI. 

AD [A0]  /dcA\  nFAD"' 
-r.)  'dr  /_  i 


(5« 


l  = 


nFAD  [A0  ]  /d?A\  nFAD1//2a1//2  nFAD1^2  [  A°  ]  a1^2  /  dc 

172  T72  V 


(M- 


^a\ 

dr  I  r. 


(591 


The  flux  relation  at  the  electrode  is  given  from  (8a) 


rdc  \  N+2 

dr")  " 

0  3  x 


rdc  \  N+2 

dr-)ro  =  jf1Al,3CB(rj,t> 


Upon  entering  the  known  boundary  conditions  and  expanding,  and 
then  substituting  these  into  (59)  we  have 

N+l 

CA(0,t)  =  Q  -  j52ajcA(rj#t) 

N+l 

cn (0 , t)  =  R  -  l  a.c_(r. 


(601 


(61; 


'B 


j=2  3  B  3't' 


with 


0  = 


R  = 


i,1/2 


A,  ,  1  nFA  [A°]D1//2n^^2  ^N+2 


*1,1 

1 


;  1/2 


IT 


A,  ,  nFA(A°)D1/2a1/2, 


*1,1 

A,  . 

Qj  = 

3  Al,l 


i 


The  discretized  equations  for  A  and  B  are  thus 


=  E  B .  -c.  (r.  t)  +  i  A.  .c Jr.  t) 

dt  _v  ij  A  j,  r-+6  .  ,  lj  A  j, 

r^-x^  1  3-1 

dcB  N+2  N+2 

dt  _  =  E  B  c  (r.  t)  +  -42L.  e  a.  .c  (r .  t) 

ri  1  j-1  B  3'  rj/h  3  =  3  13  °  D> 

The  known  boundary  conditions  from  equations  (58)  and  (61)  are 


entered, 

and  we 

have 

dcA 

N+l 

A 

dt 

r  .=X  . 

1  1 

aO^  + 

a  E 
3=2 

dCB 

dt 

N+l 

r.=X. 

1  1 

aR[  + 

a  E  : 
j=2 

where 


B.  ,Q  +  B .  -  (A.  .Q  +  A.  „_,-) 

1,1  1 ,N+2  r^+g  1,1  1 ,N+2 


2A .  .R 

b  .  ,  r  + 


i , 1  r.+g 


b..  =  -B.  .a.  +  B..  =  - (A.  .a.  -  A..) 
ID  i#l  D  iD  ri+  1,1  3  13 


The  usual  display  of  the  response  is  given  by 


E  =  EA/B"  11F  1,1  l  cA  ( 0  ,  t ) 


where  the  two  concentration  terms  in  brackets  are  found  from  equation 
(61).  The  concentration  trims  in  (01)  are  found  by  integration  of  the 
simultaneous  equation  (03)  by  the  methods  described  11,4,5,9]  previously 
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Table  1 


Time,  sec 

i  ma 

exact 

1sinulated 

1  x  10-5 

1.7104 

1.7104 

I 

o 

X 

■H 

0.5411 

0.5412 

1  X  lO-3 

0.17] 3 

0.1712 

1  x  10~2 

0.0054 

0.0054 

1  x  10*1 

0.0017 

0.0017 

1  x  10° 

0.0001 

0.0001 

ma 


Comparison  of  exact  and  -simulated  currents  for  the  simple 
change  transfer  chronoamperometric  experiment.  (N  =  6 
collocation . points . )  See  text  for  parameters. 


riGURE  CAPTION 


Spherical  electrode  geometry. 
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